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Statement of Work 


Objective 

The objective of this investigation is to evaluate the applicability of 
electromagnetic deep-sounding experituent.i using natural sources in the 
magnetosphere by incorporating Magsat data with other geophysical data. 

Approach 

The investigator shall pursue the above objective through an analysis of 
Magsat satellite data, ground-based magnetic observations, appropriate 
reference field models, and other satellite data. 

The objective will be pursued by seeking the optimal combination of 
observations which lead first to a global, and then to a regional, 
characterization of the conductivity of the Earth’s upper mantle. 

Tasks 

The following tasks shall be performed by the investigator in fulfillment 
of the above objective: 

a. Use data from Magsat satellite to constrain a long-period 
global "response function" for the average Earth at low latitudes over 
a period ranging from 6 hours to 27 days. 

b. Synchronize the Magsat data with low- latitude ground-based 
observatory data to determine the vertical gradient of the respective 
magnetic field components. Use the vertical gradient of the appropriate 
components to independently ascertain the separation of external and 
internal field contributions- 

c. Segregate the Magsat electromagnetic "response functions" 
according to the tectonic regime at the Earth’s surface and evaluate 
systematic differences between regions having lateral scale sizes on 
the order of 1000 km or greater. 

d. Theoretically evaluate problems of resolution and interpretation 
involving electromagnetic induction by temporally and spatially-varying 
magnetospher ic sources in a rotating inhomogeneous Earth as observed at 
arbitrary points in space. Use these theoretical studies to constrain 
the interpretation of Magsat data as well as to propose further applica- 
tions of satellite-based electromagnetic deep-sounding experiments, 

e. Integrate the regional response functions with other geophysical 
data in order to constrain the joint interpretation of comprehensive 
physical models. 

f. Prepare and submit to NASA periodic progress reports and a 
detailed final report documenting the results of this investigation. 


SUMMARY 


The following report describes the two stages of our present 
analysis of Magsat magnetic data which are designed to evaluate 
electromagnetic induction effects. 

The first stage consists of comparison of data from contiguous 
orbit passes over large scale geologic boundaries, such as ocean- 
land interfaces, at several levels of magnetic disturbances. The 
purpose of these comparisons is to separate induction effects 
from effects of lithospheric magnetization. The procedure for 
reducing the data in order to make such comparisons is the 
following: 1) identifying and subtracting quiet time effects; 

2) modelling and subtracting first order ring current effects; 

3) projecting an orbit track onto a map as a nearly straight line 
so it can serve as an axis on which to plot the corresponding 
orbit pass data in the context of geography. 

The second stage of our analysis consists of comparison of 
Magsat data with standard hourly observatory data. The purpose 
is to constrain the time evolution of ionospheric and 
magnetospheric current systems. Qualitative features of our 
ground based dataset are discussed. Methods for reducing the 
ground based data consist of: 1) a baseline is removed from the 
data; 2) a spherical harmonic analysis i^ performed on the data; 
3) the spherical harmonic coefficients are used to reconstruct, 
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uslnx a suitable method of interpolationy an equivalent ground 
pass” dataset corresponding to each orbit pass dataset under 
study • 


OVERVIEW 


The present analysis of Hagsat data at Brown University is 
separated into two stages. Both stages are designed to study the 
extent to which we can detect and constrain the nature of 
electromagnetic '.nduction effects on Magsat observations. In the 
first stage we examine a set of contiguous passes over specific 
land areas at specific times in a storm evolution. We utilize 
data from both minimally and highly disturbed periods as 
determined by magnetic indices and visual inspection. After 
preliminary conditioning of the data, we examine it to determine 
to what extent we can infer induction effects through correlation 
of magnetic variation with major geologic contacts. Examples o f 
such areas would be the electrical contrast across ocean*land 
interfaces; others would be the contrast between regions of high 
tectonic activity and regions which are tectonically stable. 

In the second phase of our study we compare Magsat data with 
ground based standard magnetic observatory data. We attempt to 
constrain the time evolution of ionospheric and magnetospheric 
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current eystenis before and during storm time* 

There is much motivation in the literature for searching for 
electromagnetic induction effects using satellite data (e«g*» 
Langely 1974; Didwallf 1982). Hermance (1982) points out that if 
the upper mantle is modelled as a superconductor at a depth of 
400 km» this leads to an internal to external horizontal magnetic 
field ratio of Q**34. He mentions that a storm recovery phases 
being a long period phenomenon, would be affected by the presence 
of such a conductor. He adds that for shorter periods (200sec to 
one hour), modelling the oceans as suerconductors is plausible; 
this vould lead to a value of Q»»42 at 400 km altitude. 

Carlo et. al* (1981), using ground based near-equatorial 
magnetic field data, present convincing evidence of induction 
effects from both the equatorial electrojet and Sq variations. 
The effects were noted for both quiet times and for subs torm 
events with significant enhancement of Q in the substorm data. 
He adds though that the results seemed to be highly dependent on 
the local conductivity structures of the sites involved in the 
experiment • Indeed, in an investigation in Chad and the Central 
African Republic, Farabitakoye and Mayaud (1976) asserted that the 
quiet time internal part of the electrojet field had an amplitude 
corresponding to equivalent image currents at 1200 km. 

From a theoretical model, Sugiura and Hagan (1979) show that 
the Z (vertical) component of the internal part of Sq should be 
virtually the same on the ground and at Magsat altitudes. Even 
in the dawn-dusk sector in which Magsat orbitted, they expect 
variations in Z as high as 10 gammas in some parts of the world. 
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Heraanee (1982) adds that avan tha X (horiaontal North) coaponaat 
nay be aenaitiva to ionoapharlc induction affaeta if tha 
aatellite paaaea ovar a targa acala gaologic contact. 
Eaaentially, above the current aource, the internal and axtarnal 
parta of the X component vill almoat cancel on both aidea of tha 
contact. Hovevar« the cancellation on tha lower conducting aide 
of the contact will be leaa perfect than the cancellation on tha 
more conducting aide of the contact. Thia ia bacauae of a 
difference in the equivalent homogenaoua induction paramatera on 
the two aidea of Che contact. The reault ia a amall reaidual 
field, which has, however, a pronounced characCeriaCic variation 
aaaociated with the geologic contact. For example, if one has an 
external time-varying current at 110 km with a period of three 
hours over a geologic contact with a resistivity step from 15 
ohm-m to 30 ohm-m, Hermance estimates that a horiaontal 
disturbance of 100 gammas on the ground will have a magnitude of 
about 5 gammas at 400 km, "well within the resolution of hagsat 
data . " 


II: ANALYSIS PROCEDURE 


We begin our analysis with the averaged magnetic field data 
from the Magsat Investigator-B tapes. Thia consists of data 
points at a sampling interval of about five seconds, each datum 
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being the einple average of 80 points of Chronfin fine attitude 
data. We immediately subtract main field model (4/81), supplied 
on the Investigator-*B tapes, from the averaged magnetic field 
data. We term the resulting magnetic field values "raw residual" 
data . 

IIA: SEPARATING SPURIOUS SPIKES FROM MA6SAT DATA 

We have observed short duration spike-like irregularities 
(from one to five points in duration) in the raw residual data. 
These glitches range to as large as 40 gammas in size and have 
occurred at frequencies as high as 20 incidents per orbit pass in 
passes for which data availability is supposedly complete (Langel 
et. al. 1981; p.l43, e.g. pass 1596). 

Correlation o f gl itches with bad attitude solution flags or 
shifts in instrumentation used to attain attitude data is not 
diagnostic in most cases. At times a questionable data value 
will be associated with an acceptable attitude solution while a 
"normal looking" data value directly before or after it will have 
an unreliable (motion model) or different (SCI instead of SC2 
used) attitude solution. At other times a questionable data 
value and its "normal looking" neighboring data values will all 
be associated with the same good attitude solution. In addition, 
there are numerous examples of shifts between acceptable and bad 
attitude solutions or between two different acceptable attitude 
solutions where there is no discernible accompanying spike. 

Given this complicated situation, we have decided to cull 
spurious spikes from the data by inspection of graphs of each 
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individual orbit pats. Though tha aathod ia euabaraoaa, it aaaaa 
as good as and more direct than any ayateaatic eoaputeriaed 
method involving attitude quality. It also allows one to overlay 
neighboring passes and to employ reproducibility as a criterion 
for separating attitude (or otherwise spurious) spikes from 
meaningful data. 

c 

Segments of data irregularity are reconstructed using simple 
linear interpolation from the last acceptable point before the 
irregularity to the first acceptable point after the 
irregularity. 

IIB: INTERPOLATION OF DATA TO STANDARD LATITUDE SETS 


Our objective is to compare Magsat orbit passes which are 
longitudinally coincident (or nearly so) in a quantitative sense. 
Since the latitudes at which measurements are taken are not 
precisely coincident for any two passes, we find it convenient to 
group together passes within a swath +/-2.5deg geographic 
longitude of our reference path and to then interpolate the 
magnetic data for any given pass in a particular longitudinal 
swath to a set of spatial locations at standard latitudes. The 
longitudes at these spatial locations are altered so as to nearly 
preserve the original track of each orbit pass. 


February 13th has been chosen as the quietest international 
quiet day in February. He use data from this day as a tentative 
base line when studying storm time data (see sections IV and V). 
He find it convenient to use the latitude set for each February 
13th pass as the reference latitude set for its longitudinal 
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•wath. Thit allows us to "correct" other paeaaa within the swath 
for static lithospheric and core contributions which are not 
locally accounted for in the reference field wodel (4/81) « 

Given magnetic data at latitude locations on either side of a 
standard latitude, we do a linearly weighted average between the 
magnetic data at these two latitudes (the weighting is determined 
by the relative distance of the standard latitude from each of 
the two adjacent latitudes; the longitude difference is ignored.) 
to obtain an interpolated magnetic data value at the standard 
latitude. If the magnetic value directly previous (in the sense 
of Magsat's orbit direction) to the standard latitude is padded, 
we move back a point to perform the weighted average. This 
contingency procedure may be repeated up to three times, that is, 
until we are using the third point previous to the original point 
directly behind the standard latitude. This point is somewhere 
between 120km and 160km from the point at which we want an 
interpolated value (the points are about five seconds apart and 
the satellite travels at about eight kra/sec). This corresponds 
to an equivalent earth-based separation of about 111km to 148km 


(assuming 

an altitude 

of 

500km 

). In general 

the 

maximum 

resolution 

of an 

aerial 

magnetic 

survey is limited 

to 

an 

area 

having a 

radius 

equal 

to 

or greater than about 

25Z 

of 

the 


altitude at which observations are made. This translates to a 
maximum resolution of between 81km and 138km over the lifetime 
of Nagsat. Over the course of one orbit pass in February, the 
limits are 85 km. to 125 km. Hence^ if we used data further back 
along the path to perform a weighted average, we would in effect 
be applying smoothing to our data over ranges possibly within the 
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data's resolution. Hence, if all four points previous to the 
latitude of interest have padded aagnetie values, we do not 
interpolate, but assign directly to the standard latitude, the 
aagnetic field value at the point directly after the standard 
latitude. If the magnetic data value at the point directly after 
the standard latitude is padded then in all eases the magnetic 
data value at the standard latitude will be padded. 

The above procedure is applied separately for each magnetic 
field conponent being processed. We note that one component 
might be padded in an interpolated dataset at a particular 
standard latitude while another is not if there are different 
padding patterns in the original data for the two components. 
Moreover, the two points at which data are used to interpolate e 
value for one component at a given standard latitude mey be 
different from the two points at which data are used to 
interpolate a value for a second component at the same stsndard 
latitude. This brings up the problem of how to carry along the 
attitude quality word, one of which is assigned to each data 
block (point) on the Invest igator~B tapes. We solve this problem 
by creating, for the interpolated values, a separate attitude 
quality word for each component at each point in an orbit pass. 
The attitude quality word of the interpolated component for a 
particular point is taken as the maximum (digit by digit) of the 
quality words at the two points at which the magnetic values were 
used for the interpolation of that component. 

The longitude of an interpolated magnetic data value at a 
standard latitude is obtained by taking a simple weighted average 
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of the longlt'.*dee at the two points at latitudes adjacent to the 
standard latitude in the original dataset* regardless of whether 
there are pads for any of the aagnetie data values at these two 
points. The accuracy of this Method is more than sufficient for 
preliminary investigations, especially in the mid and lower 
latitudes. For more detailed investigations and any 
investigations using data at higher latitudes where curvature in 
the longitudinal path becomes noticable, we will employ the more 
accurate method of longitude interpolation on a great circle 
path. 


Ills DEFINING QUIET TIME MAGSAT DATA 


In the context of geomagnetic field studies, a quiet period is 
most generally one during which irregular time-varying external 
field sources are at a minimum. In this case, irregular induced 
currents in the ground and the ionosphere will also be at a 
minimum (By "irregular" field sources, we refer to sources whose 
origin is associated with storm or sub-storm effects as opposed 
to sources such as the Sq current system which is a regular and 
somewhat predictable feature of the external geomagnetic field 
and seems to be unaffected by storm time activity (Regan and 
Rodrigues, 1979).). 
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Quiet orbit peteoi in Mafoet date can ba uaad to aarva two 
baaie purpoaaa. Clearly they aay ba uaad to eonatrain the 
■agnitude and location of anall aeala lithoapharie aagnatiaation 
anoaaliaa. During noiay parioda, thia avidanca can ba either 
naaked by or aiatakan for evidence of induction effacta. 
Magnetication anoaaliaa are identifiable by a repeatability in 
aeala aizct aorphology and aagnitude froa paaa to paaa« 
independent of the aaount of stora activity. So for the purpoaa 
of studying lithospheric aagnatiaation, quiet orbit pass data 
represents the signal in the data and is the subject of 
subsequent analysis. 


On 

the other 

hand, in 

the 

study 

of 

induction 

effects , 

quiet 

orbit 

pass data 

becoaes 

the 

noise 

in 

the data 

for the 

very 


reasons described in the paragraph above. In induction studies 
we are looking for transient tine-varying signals which aust be 
separated froa static nagnetic sources in the lithosphere. Henna 
for induction studies our concept of an ideal quiet day beconcs 
Bore specific than what we have outlined at the beginning of this 
section. We add the constraint that the quiet period should be 
chosen as close as possible in tine to the disturbed period we 
are choosing to study for induction effects. The reason is that, 
in addition to aagnetization anonalies, quiet pass data also aay 
contain long period non-linear trends reflecting seasonal and 
other unconstrained effects froa ionospheric/aagnetospheric 
current systeas. These effects can cause significant offsets 
between the average value of a coaponent on two different 
internatinal quiet days far apart in tine. When dealing with 
satellite data, we aust also take into account the variation of 
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altitude of dlffarant paaaea over a particular rcflott ot tha 

earth due to preeeaaion aa well aa Che deereaae in altitude over 
the life of the aisaion. Becauae of the 1/r^ dropoff of Che 

firaC decree and order terreatriallT induced parta of the field, 
there will be about a five percent difference between ita 
amplitude at 400 ka altitude and ita amplitude et 500 km 
alcicude« Thia eonatraint aakea even atronger Che ceae for uaing 
quiet paaaea cloae in atom tiae to the diaturbed period of 

intereat . 

The quiet paaa ahould alao be choaen immediately (within e 
day) before the onaet of a diaturbance being atudied. Bugiura 
(1964) mentiona that Che recovery phaae of a atora nay laat from 
several days to two weeks. It is often difficult deciding when 
Che end of a storm has merged with quiet time "noiae." This 
difficulty of deciding when the energy of the long period 

recovery phase of a storm has been expended makes clear the 
desirability of using a quiet pass immediately before, rather 
than after, a magnetic storm. 


IV: SEPARATING QUIPT TIME EFFECTS FROM MAGSAT DATA 


To study the effects of electromagnetic induction in Che solid 
earth lithospheric magnet last ion must be subtracted from the raw 
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residual magnetic field data* 

As a first approximation ve consult the Dst and Kp indices. 
Kp is a rough measure of solar corpuscular radiation activities; 
it is compiled by taking the maximum range of the most disturbed 
component over a three hour interval (after subtracting S<) 
effects) for several mid-latitude observatories. After 
correcting for known regional differences between the reference 
observatories (such as conductivity structure and seasonal 
differences), the ranges recorded at the observatories are 
referenced to a standard scale (running from zero to nine) and 
then averaged to produce one global three-hour Kp value (Regan 
and Rodriguez, 1979; Mayaud, It^SO). 

The Dst index was designed to be a measure of axially 
symmetric storm time equatorial ring current activity. It is 
compiled from the horizantal north component of the magnetic 
field at several low to mid-latitude observatories after daily 
variation effects have been subtracted out (Sugiura, 1964; Regan 
and Rodriguez, 1979; Mayaud, 1980). 

Wt. choose our initial quiet time orbit passes based on three 
criteria: low Kp, low Dst and proximity in time to magnetic 

disturbances of interest. Low Kp indicates overall low solar 
flare activity. Low Dst indicates normal ring current magnitude 
and minimal time varying ring current activity. The reasoning 
behind choosing near storm time quiet passes ie discussed in 
detail in section III. 

Based on these criteria, we have chosen the orbit passes on 
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February 13th at the quiet passes to subtract froa data of 
interest during the stora of Febraury 14th. For a given 
disturbed orbit pass d«casetf the pass on the 13th that is in the 
saae longitudinal swath (see IIB) as the disturbed data is the 
particular quiet data that is subtracted froa the disturbed data. 
The subtraction is perforaed point by point, after interpolation 
to atandard latitudes (see IIB). 


V: SEPARATING FIRST ORDER RING CURRENT EFFECTS 


After quiet tiae effects are subtracted froa a Magat orbit 
pass dataset, we should be left with only effects froa stora 
linked tiae-varying external current systeas and their induced 
terrestrial counterparts. The aost proainent of these effects is 
the enhanceaent of the east-west ring current systea. We wish to 
reaove a first order fit of the ring current effect in our 
aagnetic data in order to 1: study the second order aorphology of 
the stora tiae enhanced ring current; 2; study effects froa 
localized finite-source systea:. unaasked by the reaoval of the 
first order ring current field. 

We aodel the ring current systea with a siaple two paraaeter 
function ; " A sin (0) + B cos (9) ", "6" is geographic latitude 

here. This allows for a siaple sinusoid with an asyaaetry with 
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respect to the geoRrephic equator. To avoid contaalnatlon of the 
fit by auroral currents we only perfora the fit on data for 
-45deg < 3 <45deg. The resulting model however is subtracted froa 
the complete half-pass (ascending(dusk) or descending(dawn) 
node). We compute separate models for the daawn and dusk 
segments of each orbit pass. This is done in order to identify 
any first order asymmetries between the dawn and dusk portions of 
the ring current either in amplitude or latitude position. 

Properly, a first order latitude asymmetry, specifically a 
tilt in the ring current about its equatorial axis, should be 
modeled by a function With an '*exp(i* 6 )" longitude dependence. 
However, as our data is divided naturally, pass by pass, into 
areas of small longitudinal extent, we ignore this dependence in 
our preliminary studies. Also, modelling the ring current is 
more properly done using geomagnetic coordinates. As we allow 
for latitude assymetry of the sinusoid however, this becomes 
academic (in stage two of our analysis we do rotate all our 
coordinates to refer to the geomagnetic pole). 

There are approximately 280 points in a quarter orbit pass 
(~45deg < 9 <45deg).' Before we perform our sinusoidal fit, we 
cull the data by means of its attitude quality. Our requirement 
is that 1: a Quest solution (Langel et.al., 1981) for the 
attitude must have been computed (attitude flag digit d <* 3); 2: 
the Quest solution residual be less than 20 arcseconds (attitude 
flag digit b < > 3). If this leaves us with too sparse a dataset 
for a particular fit, we relax the requirement for a Quest 
solution at every point and allow data points with a motion model 
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attitude eolution for up to three date pointe in a row following 
a point with a Quest attitude solution with an acceptable 
residual. This is well within the half~ainute liait for use of 
action aodel results recoaaended by NASA (Langel et.al., 1981; 
p.3). However, so far, we have found this unnecessary, as 
anywhere between 80 and 100 points have been available for the 
fit for each orbit pass dataset. 


VI: VIEWING ORBIT PASS DATA IN CONTEXT OF GEOGRAPHY 


Once our selected passes of Magsat orbit data are processed as 
detailed in sections I-V, they are ready for examination for the 
effects of electromaanetic induction within specific large scale 
conductivity anomalies within the earth. In order to make a 
qualitative assessment at satellite altitudes of induction 
effects associated with specific geologic areas, we present an 
orbit pass dataset against a backdrop of the geography underlying 
the orbit pass. We then examine correlations between edges of 
large scale conduction anomalies (such as the ocean-*land 
interface) and satellite magnetic field variations. 


To do this we employ a computer mapping program able to plot a 
map of the world, from a digitized map database, using the 


Mercator projection with zero degrees true scale latitude (Hinks, 
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1921). Upon this map we project an orbit patht given a space 
series of angular positions that represent the orbit path. A 
Magsat orbit path however* being near-polar in trajectory* 
appears as a bell-shaped curve in this type of projection. It 
would thus be an inconvenient axis over which to plot a 
corresponding orbit pass dataset. 

To improve this situation we have developed computer code 
modifying the mapping program to allow an arbitrary rotation of 
the projection. The rotation is specified by choosing a 
particular great circle path as the great circle of true scale 
for the projection. Upon fitting a great circle to an orbit 
path* we choose that great circle as the great circle of true 
scale determining the projection rotation. The resulting 
projection depicts the orbit path as a nearly straight line. It 
is not completely straight because* as a result of the rotation 
of the earth beneath the satellite, an orbit path is not a great 
circle in earth-fixed coordinates. Figure 1 shows five 
neighboring orbit passes in a longitudinal swath with a 
projection rotated beneath the orbit pass in the center of the 
swath. The continental boundaries are easily recognizable in the 
vicinity of the orbit passes. 

Thus our two objectives are met; 1) An orbit path is projected 
as a nearly straight line so as to serve as an axis on which to 
plot the corresponding orbit pass data; 2) The surrounding 
geography is plotted in a relatively non-distorted manner so that 
correlations between magnetic variations and large scale 
conductivity anomalies can be made. Figure 2 shows Magsat 
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■tgnetlc field date from one orbit pace proeeesed at described in 
sections I-V. The horiaontal scale length is the ease as the 
horizontal scale length of the orbit pass tracks plotted in 
Figure 1 . 


VII: PREPARING GROUND BASED DATA AND NAGS AT DATA 
FOR INTERCOMPARISON 

VI lA: GROUND BASED OBSERVATORY DATA 

Our ground based data base consists of hourly means for the 
Magsat epoch from about forty observatories (we are expecting 
data from twenty to forty additional sites to be added to our 
database from WDC-A within the next few months.). In our initial 
correlations between ground and satellite data, we will 
concentrate on the month of February. Figure 3 shows data 
plotted from three representative observatories at widely spaced 
locations around the earth for this time period. One notes 
immediately the wide variety of well-defined phenomena available 
during this month. In particular, the magnetic storm on February 
14th has a clear main phase with a classic recovery phase of some 
five days in extent. As an aside we note the problem we would 
have choosing February 21 or 22 as a quiet day baseline. 
Although both are among the quietest days in the month, they are 
unfortunately superimposed upon the exponentially decaying 
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recovery phase of the February 14 stora and would significantly 
bias any long period inforaation we aight have in our data. 

The stora of February 6 also has ‘possibilities for analysis. 
It is convenient that we have two storas of respectable aaplitude 
so close together and hence with roughly a coaaon baseline. 

An interesting phenoaenon is apparent froa just a preliminary 
visual inspection of the data. At low latitudes and particularly 
in the southern hemisphere in general, the main phase of the 
February 14 storm is significantly the stronger of the two 
storms. Specifically. 802 of the observations at geomagnetic 
latitude 45deg or below, and all observatories in the southern 
hemisphere (ranging to -43deg geomagnetic latitude) have 
significantly stronger February 14 storm amplitudes. On the 
other hand, 532 of the observatories located above 45deg 
geomagnetic latitude have a main phase amplitude for the February 
6 storm that is greater or equal to the main phase amplitude of 
the Febrary 14 storm. 


There are also strong local time correlations. All of the 
sites below 45deg geomagnetic latitude that had stronger February 
6 storm amplitudes lie in a single 80deg sector between 55deg and 
135deg geomagnetic longitude. None of the remaining sites below 
45deg geomagnetic latitude fall into this sector. It would 
obviously be of some use to correlate these patterns with those 
constraints on ionospheric current flow that we deduce from our 
comparison of earth and satellite data. 

The February ground database also has a nice set of quiet days 
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suitable for Sq study. The five day span, February 1~5* contains 
three of the 10 quietest days of the aonth (Langel et. al.y 
1981). 

Before we can compare ground based data to Nagsat raw residual 
data, we must remove a baseline from the ground based data. 
Subtracting the annual mean from the data at each sits is not 
sufficient; it ignores any long period, such as seasonal, 
variations; these are corrected for in the Magsat 4/81 spatial- 
temporal main field model supplied on the Invest igator-B tapes 
(Magsat Information Bulletin, 1981). We employ a method whereby 
we compute the baseline at each site by taking a simple average 
of the three hourly values about local midnight from each of 
three quiet days scattered through the month that appear 
uncontaminated by recovery phase effects. Though not as refined 
as the continuous temporal correction employed in Magsat (4/81), 
our method does seek to minimize seasonal and recovery phase 
effects. 

Our reduction of the ground based hourly data continues with 
an hour by hour spherical harmonic analysis (SHA) of the data in 
the time domain. The SHA is taken up to degree and order three 
and uses only those sites that fall into the geomagnetic latitude 
range -45 to <*-45 degrees. The power at auroral latitudes at 
spatial wave numbers higher than degree three is strong enough to 
cause significant aliasing to any low order SHA fit we perform on 
time domain data. The auroral and polar cap observatories are 
therefore omitted from our initial analysis. 
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The SNA will leave ua with a set of coefficients at each hour 
in February froai which to reconstruct the ground aagnetic field. 
Our goal is to reconstruct an equivalent space series underneath 
the great circle orbit path of a given pass using these spatial 
harmonic coefficients of the ground based data. Obviously we 
will not be able to represent any feature of the anomalous ground 
magnetic field with a spatial wave number n>3. Because of this, 
and because the SHA was not constrained by data at high latitudes 
in any case, we will only attempt a reconstruction of the field 
under an orbit pass for the geomagnetic latitude zone between >44 
and <^45 degrees. 

For a given Magsat orbit pass, we generate a ground magnetic 
value underneath every Magsat data point. If the value at that 
point was measured right on the half hour (tabulated hourly 
values are evaluated every hour on the half hour), say 09:30, we 
use directly the coefficients from the SHA at 09:30 to 
reconstruct the magnetic field at the latitude and longitude of 
the Magsat data point. To generate a ground data value 
corresponding to a Magsat data value not measured precisely on 
the half hour, we perform a bi'cubic spline interpolation using 
the three nearest neighbor sets of coefficients from the SHA. 
Using this method we generate an equivalent “ground pass" dataset 
for each orbit pass dataset in February on the geomagnetic 
latitude interval from -45 to +45 degrees. 
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FIGURES 


Fixure 1 : Orbit passes in the same lonxitudinal swath from 

February 13, 15, 19 and 21. The tracks and the map are plotted 

using a rotated Mercator projection so that the rotated aero 
latitude and zero longitude point of the map is the point of 
northernmost approach of the February 17 pass. This tends to 

minimize the curvature of the satellite trajectory relative to 
conventional Mercator projections. 


Figure 2 : Data from Magsat orbit pass #1627, February IS. 
Only data between -45 degrees and 45 degrees geographic latitude 
is shown. The data has had attitude spikes, quiet time effects 
and first order ring current effects removed. 


Figure 3 : Ground based standard observatory hourly values 
from Novosibirsk (45deg, 158deg), Memambetsu (34deg, 209deg) and 
Port Moresby (-19deg, 218deg). The coordinates are geomagnetic 
latitude and longitude. 
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TIME UNITS (MIN) 

in.n 2S,l 31.3 37.6 «3.9 SO. 1 M.« K.7 M.O 75.3 
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